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. We are responding to the comment of Tyshetskiy and Vladimirov [lj (hereafter referred to as TV) on a recently 
published paper by Shukla and Eliasson (hereafter referred to as SE), with the corresponding Errata regarding 
the discovery of a novel attractive force on stationary test ion charges of the same sign shielded by degenerate electrons 
in a dense quantum plasma. For their purposes, Shukla and Eliasson [2[ used the lineraized quantum hydrodynamic 
(LQHD) equations composed of the electron continuity equation, the inertialess non-relativistic electron momentum 
equation with electrostatic and quantum forces [e.g. due to the quantum statistical pressure that accounts for electron 
i— | , degeneracy, the quantum recoil effect associated with electron tunneling through the Bohm potential via overlapping 
tIh ' of quantum electron wave functions and weak- localization of degenerate electrons, and electron-exchange and electron 
SpH, correlation (EXC and EC) effects due to electron its spin together with the Poisson's equations. In Fourier space 
the latter yielded the SE potential distribution (Eq. 2 in TV) due to the dielectric constant that resulted by the 
quantum forces that we described above. TV in their recent evaluation have come up with a misnomer that the 
short-range (of atomic scales) SE potential is insignificant effect compared to the long-range (in comparison with the 
Fermi-Thomas length Aft = fFe/(v3a>p e ), where VFe — hkF/m e is the nonrelativistic Fermi electron speed, h, the 
Planck constant divided by 2%, &f = (S^rio) 1 ^ 3 the Fermi wave-number, uq the unperturbed electron number-density, 
ujpe = (47rnoe 2 /m e )'l/2) the electron plasma frequency, e the magnitude of the electron charge, and m e the rest mass 
of electrons) Friedel oscillations (FOs), predicted by the Lindhard's dielectric constant e(0, k) [J] in the zero-frequency 
and zero-Fermi electron temperature limits, by ignoring electron-exchange and electron correlation effects, and which 
' appears as Eq. (4) in TV, where k is a non-normalized wave number. The FOs are caused by discontinuity in the 
Lindhard'd dielectric function at the Fermi-surface, where k = 2kp. In their Fig. 1, TV erroneously compare the 
SE potential profile deduced from their Eq. (2), depicted by a dashed line, which includes electron-exchange and 
t-H \ electron correlation effects, with their oscillatory potential profile (depicted by the solid line) that is obtained by the 
^ 1 numerical integration of their Eq. (3) which uses Eq. (4) that completely neglects EXC and EC effects. Such a 
comparison in Fig. 1 of TV is meaningless because the SE potential distribution (SEPD) profile, as depicted in Ref. 
[H significantly differs from the dashed profile in Fig. 1 of TV, besides the fact that the SEPD profile encompasses 
\ the physics of the EXC and EC effects that are not included in the Lindhardt's dielectric constant (Eq. 4) which has 
been used In Eq. (3) of TV. 

In the following, we present a corrected analysis that TV gave in their earlier comment [l[ , and demonstrate that the 
SE potential distribution without the EXC and EC effects (referred to as the reduced SEPD) is a different phenomenon 
which has nothing to do with the FOs. It should be stressed that both the generalized SEPD and the present reduced 
SEPD distributions do not assume /cAf *C 1 or k = 2kp. For our purposes, we shall use the LQHD dielectric constant 
• \ from Ref. [2j without the EXC and EC effects, which is of the form 
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and the dielectric constant (in the zero- frequency limit) of Lindhard Q , which uses the random- phase-approximation 
(RPA) for noninteracting completely degenerate electron fluids, and which reads [7| 
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where w± = ±frk 2 /2m e . It is interesting to note that Eq. (1) is exactly the same as the one obtained from the 
quantum kinetic theory's dielectric constant [1] [viz. the reduced Eq. (14)] in the zero-frequency and zero Fermi 
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electron temperature limits, without any approximation on k\p. Consequently, the electrostatic potential around a 
stationary test ion charge, in dimensional from, reads as 
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where 3 stands for the imaginary part of integral. It should be noted that the dielectric function can be either 
D® HD (0, k) or D RPA (0, k) given by Eqs. ([IJ and ©, corresponding to the SE or TV schemes, respectively. 

It is well known Q that the Lindhard dielectric constant, given by Eq. (2), gives rise to the Kohn anomaly Q due 
to a singularity at the specific wavelength Xk = ir/kp at the electron Fermi surface. Such a singularity is known to 
produce the periodic local density pattern of the form n(r) ~ sm(2kpr)/r d at distance r around a stationary test 
charge [l(| far from the charge itself. This is referred to as FOs in d-dimensional compact Fermi electron-sea. In 
their coment, TV [H have presented their numerical simulation results based on the RPA dielectric function (Eq. [2]), 
claiming that the amplitude of SE potential distribution (shown by the solid line in Fig. 1 of TV) is comparable 
or even less than that of FOs. Our objective here is to point out an obvious mismatch in TV's calculation and 
inappropriate comparison between the reduced SEPD profile and the potential profile deduced from their Eq. (3) by 
using Eqs. (4) and (5) which lack essential interaction electron picture. In order to evaluate the outcome of the two 
approaches on equal footings and to be consistent with TV's scaling, we use the same normalization as in Ref. 
Thus, fc _1 and r are in units of Xp. Accordingly, one can cast Eqs. (p} and ([5} for both the reduced SE (denoted by 
QHD) and RPA model which ignore the EXC and EC effects. We have 
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where a = (3hu) pe /2m e v Fe ) 2 = (9/7r 5 ) 1 / 3 (r /a_B) = (9/n 5 ) 1 ^ 3 r s (r s is the Brueckner parameter), for which r - 
and as = h 2 /e 2 m e are the inter-electron spacing and the Bohr radius of an hydrogen atom, respectively. We note 
that a here is different from the a in Ref. [2J because of the absence of the electron exchange and electron-correlation 
potential (V xc ) in the present analysis. 

The normalize potential, (j>(R), defined based on Eq. ([3]) in the TV-compatible scaling, may be written as 



where the electrostatic potential <fi(R) is normalized by Q/Xf, with Q being the ion charge state, and D(0,k) the 
static dielectric constant for the reduced SE and RPA cases. 

Figure I depicts the inverse of the dielectric constant, 1/D(0, k), for the reduced QHD and RPA cases for a given 
electron number density no = 9.63 x I0 22 cm~ 3 , which corresponds to a ~ 2.06 (r s ~ 6.68) used in Fig. I of TV's 
comment [l|. It is noted that in the long- wavelength limit 1/13(0, oo) = 1 is fulfilled for both cases, as expected. 
The Kohn-anomaly due to the singularity in the Fermi electron- liquid dispersion embedded in Eq. (2), which leads 
to weak density oscillations (FOs) around an impurity test charge, occurs exactly in the middle region of the dashed 
curve (the RPA case) at k = 2kp ei where it seems to be broken. 

Figure 2, on the other hand, exhibits the variation of the potential profiles for the reduced QHD (the left plot) and 
RPA (the right plot) cases for different electron number densities of our concern in this reply. These profiles from 
thick to thin curves (the highest to the lowest electron number densities) correspond to a ~ {0.19, 0.93, 2.06, 2.5} or 
r s ~ {0.62,3,6.68,8.1}, respectively. The potential distribution profile corresponding to the electron number density 
of no = 9.93 x 10 22 cm -3 (denoted by an asterisk symbol in the plots) is the one given by TV. Indeed, the expected 
oscillations from the RPA theory are evident from Fig. 2. However, while there appears to be a pattern for the 
electric potential depth and its location in the reduced QHD case, in the RPA case no such pattern is found at 
least for characterizing the first minimum in the oscillations. Such lack of consistent pattern can be marked as a 
disadvantage for a theory. While TV fail to simulate the reduced QHD potential minimum location correctly in their 
paper, they also underestimate its depth as is apparent when comparing Fig. 2 for no = 9.93 x 10 22 cm~ 3 case here 
with that of Fig. 1 in TV's comment paper [lj. The first positive local minima of the corresponding, a — 2.06, the 
potential profile for the RPA case is not even shown in Ref. [1] . The location of this positive minimum (not shown in 
Ref. [l|) is around the minimum of the corrected potential profile for the reduced QHD case (e.g. see Fig. 2 here) for 
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FIG. 1: (Color online) The inverse of the dielectric constant versus the parameter K using the reduced QHD (solid) and RPA 
(dashed) cases for a = 2.06 (n = 9.63 x 10 22 cm" 3 ) corresponding to Fig. 1 of TV. 
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FIG. 2: (Color online) The profile of the potentials as a function of R for different values of the electron number density no for 
the reduced SE (QHD model) (left) and RPA (right) models for different electron number-densities no = (7.94, 9.63*, 21.3, 104) x 
10 22 cm -3 (the star mark denotes the number-density used in TV's comment) shown by (Brown, Red, Blue, Black) colors, 
respectively. Similar colors in the two models correspond to the same electron number density. 



the same electron number density. It is found that for all electron number density values used here (including the one 
used Ref. [l[) the amplitudes of the RPA profiles are lower than that of the reduced QHD case, which is in contrast 
to the findings of Ref. [l| . Only a somewhat qualitative agreement between the reduced QHD and RPA cases seems 
to appear in Fig. 2 for the lowest metallic density of no = 7.94 x 10 22 cm~ 3 and the highest electron number density 
no = 104 x 10 22 cm -3 , with the Thomas-Fermi-like charge screening. 

On the other hand, it has to be mentioned that the results of the RPA case without considering the important 
electron interaction features, such as EXC and EC effects, are highly unreliable. TV claim that such an extension 
of the RPA theory is being elaborated by them. However, such an extension has in fact been already carried out by 
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Singwi et.al. [ll, 12 1, who have proven that the Lindhard dielectric constant 0,0]; obtained by using the RPA, is 
only adequate for the description of the plasmon excitation modes in the long-wavelength limit and its applicability 
is limited to the high electron number density. This theory has also been shown to lead to nonphysical negative 
pair-distribution function for the entire range of the metallic density, which is caused by the lack of the local field 
correction and short-range effects in the RPA theory [ll[ . The screening density (which is proportional to the screening 
electrostatic potential) for the RPA and improved Hubbard and Singwi models has been shown for r s — 3, 6 values 
in Figs. 6 and 7 of Ref. [ljj]- At a first glance of these Figs, it appears that the improved screening tends to deepen 
the first valley making it look more like those of the reduced QHD rather than those of the RPA, as shown in Fig. 



2. Let us closely inspect Figs. 6 and 7 of Ref. [l2j against the potential profiles displayed in Fig. 2 for the electron 



number densities no — 9.63 x 10 22 cm~ 3 and no = 21.3 x 10 22 cm~ 3 corresponding to Brueckner parameters of r s = 3 
and r s = 6.68, respectively. It is observed that in improved models while the location of the principal minimum is 
almost unchanged, its depth increases with a decrease in the electron number density. These findings are in complete 
agreement with the reduced QHD findings for the electric potential profile, as shown here in Fig. 2. It also sharply 
contrasts the findings of the RPA result that is presented in our Fig. 2. It is evident that the noninteracting RPA 
model used in Ref. [l| does not lead to physically consistent results. For metallic densities of no = 9.63 x 10 22 
cnr 3 and no = 21.3 x 10 22 cm -3 it does not even present a well pronounced principal negative potential valley, a 
characteristics of the quantum plasma screening. Such a well pronounced negative potential valley can even be found 
in density functional theory (DFT) simulation results [l3| of quantum plasma systems, which deepens the potential 
well by including local density corrections. 

In summary, we conclude that only a quantum kinetic theory which includes all quantum forces, as described in 
Ref. must be used to obtain physically consistent results for the quantum screening. The RPA model used 
by TV is not consistent with the LQHD, DFT and improved Hubbard and Singwi models and should therefore be 
discarded. There are oscillations in the screened potential tail of a quantum electron fluid due to the Kohn-anomaly, 
insignificant compared to the well-pronounced reduced SE effect. Such insignificant oscillations can also be probed 
by extra zooming-in of the reduced SE potential's tail, which might appropriately account for the FOs phenomenon. 



[1] Yu. Tyshetskiy and S. V. Vladimirov, Comment on "Novel Attractive Force between Ions in Quantum Plasmas", submitted 

to Phys. Rev. Lett, (under review); see also E-print: arXiv: 1212.4286 ^1 [physics. plasm-ph] 19 Dec 2012. 
[2] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 165007 (2012). 

[3] P. K. Shukla and B. Eliasson, Phys. Rev. Lett. 108, 219902 (E) (2012); ibid. 109, 019901 (E) (2012). 
[4] J. Lindhard, Kg. Danske Videnskab. Selskab, Mat. Fys. Medd. 28, No.8, 1 (1954). 

[5] Yu. L. Klimontovich and V. P. Silin, Dokl. Akad. Nauk SSSR 82, 361 (1952); Zh. Eksp. Teor. Fiz. 23, 151 (1952). 
[6] D. Bohm and D. Pines, Phys. Rev. 92, 609 (1953). 

[7] E. M. Lifshitz and L. P. Pitaevskii, Physical Kinetics (Butterworth-Heinemann, Oxford, 1981). 

[8] N. L. Tsintsadze and L. N. Tsintsadze, Europhys. Lett. 88, 35001 (2009). 

[9] W. Kohn, Phys. Rev. Lett 2, 393 (1959). 
[10] George E. Simion and Gabriele F. Giuliani, Phys. Rev. B 72, 045127 (2005). 
[11] K. S. Singwi, M. P. Tosi and R. H. Land, Phys. Rev. 176, 589 (1968). 
[12] K. S. Singwi, M. P. Tosi and R. H. Land, Phys. Rev. B 1, 1044 (1970). 
[13] Jeng-Da Chai and John D. Weeks, J. Phys. Chem. B 108, 6870 2004. 



